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Abstract 

Background: The shape of grass leaves possesses great value in both agronomy and developnnental biology 
research. Leaf rolling is one of the innportant traits in rice {Oryza sativa L) breeding. MYB transcription factors are 
one of the largest gene families and have important roles in plant development, metabolism and stress responses. 
However, little is known about their functions in rice. 

Results: In this study, we report the functional characterization of a rice gene, OsMYBIOSL, which encodes an 
R2R3-1\/1YB transcription factor. OsMYBIOSL was localized in the nucleus with transactivation activity. Overexpression 
of OsMYBIOSL in rice resulted in a rolled leaf phenotype. Further analyses showed that expression levels of several 
cellulose synthase genes {CESAs) were significantly increased, as was the cellulose content in OsMYBIOSL overexpressing 
lines. Knockdown of OsMYBIOSL by RNA interference led to a decreased level of cellulose content and reduced 
mechanical strength in leaves. Meanwhile, the expression levels of several CESA genes were decreased in these 
knockdown lines. 

Conclusions: These findings suggest that OsMYBIOSL may target CESA genes for regulation of cellulose synthesis 
and could potentially be engineered for desirable leaf shape and mechanical strength in rice. 

Keywords: OsMYBIOSL, Leaf rolling, MYB transcription factor. Cellulose 



Background 

Leaves are the main organs of photosynthesis and play a 
crucial role in plant development. Leaf rolling is a fre- 
quently observed phenotype in fields and rice mutant 
populations [1-3]. The appropriate leaf rolling has been 
regarded as one of the most important elements in rice 
ideotype breeding as it improves photosynthetic efficiency, 
accelerates dry-matter accumulation, reduces solar radi- 
ation on leaves, decreases leaf transpiration under drought 
stress and raises grain yield [2,4-9]. However, severe leaf 
rolling could lead to growth retardation, developmental 
defects and crop yield reduction. 
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To date, many studies have been performed to 
characterize the genes controlling leaf rolling by analyses 
of the rolled leaf mutants in rice. Twelve rolled leaf (rl) 
mutants in rice are reported and the corresponding 
genes are mapped genetically, of which, six rl genes are 
mapped on rice chromosomes by morphological markers 
and the rest are directly mapped in rice genome by mo- 
lecular markers [10-14]. Among these 12 rl mutants, r/9 is 
the first to be cloned and analyzed in detail. ROLLED 
LEAF9 (RL9)ISHALL0T-LIKE1 {SLLl) encodes a KANADI 
transcription factor, and rl9/slll mutants display extremely 
incurved leaves due to the defective development of scler- 
enchymatous cells on the abaxial side [1,15]. Some other 
genes are also found to be related to leaf rolling in rice. 
Loss-of-function of OsCSLD4 gene, which encodes a 
cellulose synthase-like protein, results in phenotypes of 
reduced leaf width and semi-rolled leaves, possibly due 
to the significantly smaller bulliform cells in mutants 
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[16-19]. Loss-of-function of CONSTITUTIVELY WILTEDl 
{OsCOWl) I NARROW LEAF? {HALT), encoding a mem- 
ber of the YUCCA protein, causes rice leaves to roll up- 
ward [20,21]. Loss-of-function of CURLY FLAG LEAF 
(CFLl), encoding a WW domain protein, induces rolled 
leaves which have a wrinkled epidermis and contain pa- 
pillae with reduced wax content [22]. ROLLED LEAF 14 
{RL14), encoding a 20G-Fe (II) oxygenase family pro- 
tein, is reported to modulate rice leaf rolling by affecting 
secondary cell wall formation and buUiform cell devel- 
opment [2]. Mutation of SEMI-ROLLED LEAF 1 (SRLl), 
encoding a putative GPI-anchored protein, results in 
leaves rolled adaxially in rice due to the increased num- 
ber of buUiform cells at the adaxial cell layers [23]. In 
addition, recent studies show that overexpression of rice 
gene OsAG07, which encodes an Argonaute (AGO) 
family member, results in the leaf blades curling upward 
[3]. OsHBl is a member of the Class III homeodomain 
leucine zipper family of genes, overexpression of its micro- 
T^A/Ai 66- resistant version could result in adaxially rolled 
leaves due to the reduced sclerenchyma and formation of 
buUiform ceUs on the abaxial side [24] . Taken together, the 
major factors possibly affecting leaf rolling include the al- 
teration of osmotic pressure and/or turgidity in buUiform 
ceUs, the impairment of ceU polarity establishment and 
ceU differentiation. However, the mechanisms underlying 
leaf rolling in rice and, particularly, the involved transcrip- 
tion factors remain largely unknown [1]. 

MYB proteins belong to a super-family of transcrip- 
tion factors, which prevail in all eukaryotes and share a 
highly conserved DNA-binding domain, MYB [25-27]. 
Depending on the number of adjacent repeats in the 
MYB domain, MYB proteins can be classified into four 
subfamiUes, including IR-MYB, R2R3-MYB, 3R-MYB, 
and 4R-MYB [26,28]. Among them, R2R3-MYB subfam- 
ily is a typical representative one, whose members con- 
tain two MYB repeats in their DNA-binding domains 
[29,30]. MYB proteins in plants are involved in a variety 
of plant-specific processes, including primary and secondary 
metabolism, cell fate and identity, developmental processes, 
and responses to various stresses [25-27,31-40]. Genome- 
wide analysis shows that there are approximately 183 MYB 
genes in rice, of which 109 encode R2R3-MYB proteins 
[41,42]. So far, only a few number of MYB proteins, includ- 
ing OsGAMYB, OsMYB2, OsMyb4, OsMYB3R-2, CSA, 
MYBS3 and OsMYB2P-l [43-49], have been functionally 
characterized. The functions of most of MYB proteins are 
stUl unknown in rice. 

To characterize functions of MYB transcription factors 
in rice, we overexpressed several MYB genes in Kasalath, 
an indica cultivar, with the transgene constructs contain- 
ing the full-length cDNAs of rice MYB genes, driven by 
maize ubiquitin promoter. Of the transgenic lines, one 
line overexpressing the full-length cDNA of Os08g05520, 



designated as OsMYBlOSL hereafter, displays a rolled leaf 
phenotype. OsMYBlOSL encodes an R2R3-MYB transcrip- 
tion factor. Our study shows that it localizes in the nucleus 
and possesses transcriptional activity. We detaU the 
phenotypes of OsMYBlOSL overexpressing (OE) and RNA 
interference (RNAi) knockdown plants, including the al- 
tered leaf shape, the changed cellulose content, and the im- 
paired mechanical strength. The roles of OsMYBlOSL in 
leaf shape formation and ceUulose synthesis are discussed. 
We propose the potential application of OsMYBlOSL in 
molecular breeding of rice. 

Results 

OsMYBlOSL encodes an R2R3-MYB transcription factor 

To discover transcription factors controUing leaf devel- 
opment, we screened the rice lines ectopically expressing 
rice MYB genes under the control of maize ubiquitin 
promoter. One line overexpressing LOC_Os08g05S20 
was selected for further study due to its particular leaf 
shape, such as upward curling of the leaf blade. Accord- 
ing to the rice genome annotation database (http://rice. 
plantbiology.msu.edu), LOC_Os08g05520 encodes a puta- 
tive R2R3-MYB family transcription factor with a length 
of 359 amino acids and a molecular mass of approximately 
40 kD. The Pfam database (http://pfam.sanger.ac.uk/) 
shows that the deduced protein has two MYB DNA- 
binding domains (PF00249) at the N-terminus (Figure lA). 
As revealed by phylogenetic analysis of the related MYB 
transcription factors in Arabidopsis thaliana and rice, 
Os08g05520 is closely related to Atlg63910 (AtMYB103) 
[50] (Figure IB). Protein sequence alignment showed 
that they are highly conserved in the predicted R2- and 
R3-MYB DNA-binding domains (Figure lA). We hereby 
designated Os08g05520 as OsMYBlOSL {Oryza sativa 
MYBlOSLike), 

Nuclear localization is one of the significant features 
of transcription factors. To determine its subceUular 
locaUzation, we fused green fluorescent protein (GFP) 
to the C-terminus of OsMYB103L to produce an 
OsMYB103L-GFP fusion protein. We monitored the 
fluorescence of the transiently expressed fusion protein 
in rice protoplasts and onion (Allium cepa L.) epider- 
mis cells. In both cases, the fluorescence signals of the 
fusion protein were observed predominantly in nuclei 
(Figure 2 A and Additional fUe 1: Figure SI). WhUe in 
the GFP alone control, fluorescence signals were observed 
in nuclei and cytoplasm (Figure 2 A and Additional file 1: 
Figure SI). These results indicate that OsMYB103L is a 
nuclear-localized protein. 

Transcriptional activation is a characteristic feature of 
many transcription factors. The full-length, N-terminal 
MYB DNA-binding domain and C-terminal putative 
activation domain were fused in-frame with the GAL4 
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Figure 1 OsMYBIOSL (Os08g05520) is an R2R3-IVIYB transcription factor and iias tlie most sequence iiomology to AtlVIYB103 
(At1g63910) in Arabidopsis thaliana. (A) Protein sequences alignment of OsMYBlOSL and AtMYBlOS. OsMYBlOSL and AtMYBlOS are liighly 
conserved in predicted R2- (blacl< line) and R3-MYB DNA-binding domains (grey line). (B) Phylogenetic relationships among OsMYBlOSL and related 
R2R3-MYB proteins in rice and Arabidopsis tlialiana. An unrooted phylogenetic tree was generated with the full-length amino acid sequences. 



DNA binding domain, respectively. We conducted a 
transcriptional activation analysis in yeast by using the 
fusion proteins. The analysis showed that the full- 
length and C-terminal region of OsMYBlOSL activated 
the expression oi AURl-C and MELl reporter genes in 
yeast (Figure 2B), indicating that OsMYBlOSL has the 
capability of transcriptional activation. 



Expression of OsMYB103L varies in different organs 

We determined the expression level of OsMYBlOSL in 
various rice organs, including root, culm, leave, and 
panicle using quantitative real-time PGR (qRT-PCR). 
The expression level of OsMYBlOSL was highest in 
culm, lower in root, leaf and mature panicle, and lowest 
in young panicle (Figure SA). 
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Figure 2 Subcellular localization and transactivation analysis of OsMYB103L (A) Subcellular localization of OsMYBlOSL GFP and 
OsMYB103L-GFP fusion gene under the control of the CaMV35S promoter were expressed transiently in rice protoplasts. Left to right: GFP 
fluorescence image, transmission image and merged image. Bar = 2 |jm. (B) Transactivation analysis of different regions of OsMYBlOSL fused with 
the GAL4 DNA binding domain in yeast. The full-length, N-terminal MYB DNA-binding domain (1-160 amino acids) or the C-terminal putative 
activation domain (161-359 amino acids) were respectively cloned into pBD-GAL4 vector containing AURl-C and MELl reporter genes and then 
transformed into yeast host strain Y2HGold. The pBD vector was used as a negative control. The expression of AURl-C confers strong resistance 
to the highly toxic drug Aureobasidin A (Aba). The activities of a-galactosidase encoded by MEL-1 were examined by X-a-gal staining. 



In addition, the expression pattern of OsMYBlOSL was 
examined by the ^-glucuronidase (GUS) reporter gene 
driven by OsMYBlOSL putative promoter in rice plants. 
Examination of GUS activity in three independent trans- 
genic lines revealed that the GUS signals were stronger 
in culm and weaker in vascular bundles of coleoptile, root 
and leaf (Figure 3B to 31). GUS staining was stronger in 



mature panicle than young panicle (Figure 3F and 3G). 
Further examination showed that GUS activity was detect- 
able in pedicel, lemma and palea while no GUS activity 
was detected in sterile lemma, lodicule, pistil, and stamen 
in mature panicle (Figure 3H and 31). These results indi- 
cate OsMYBlOSL has diverse expressions in various tissues 
and organs in rice. 
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Figure 3 Expression analysis of OsMYBWSL (A) Expression of OsMYBlOSL in various organs determined by qRT-PCR analysis. Roots were 
liarvested from 2-weel<-old plants. Culms, leaf blades and mature panicles were harvested from rice plants at heading stage. Young panicles were 
collected at the 5 cm stage. Rice ACTINl gene was used as internal control. Error bars represent the SD of transcript levels determined from three 
independent replicates. (B) to (I) Examination of GUS activity in the transgenic plants expressing OsMYBlOSL pro::GUS. The GUS activity is shown 
in coleoptiles (B), root (C), culm (D), leaf (E), young panicle (F) and mature panicle (G-l). GUS activity was also detected in pedicel, lemma and 
palea in mature panicle (H), while no GUS activity in sterile lemma, lodicule, pistil and stamen (I). 



Overexpression of OsMYB103L results in a leaf rolling 
phenotype 

Twelve independent transgenic lines overexpressing 
OsMYBlOSL (OE) were obtained. The expression level 
of OsMYBlOSL was determined by qRT-PCR (Figure 4A). 
Three independent transgenic lines (OE-1, OE-4 and 
OE-8) with the highest expression levels were selected 
for further study (Figure 4A). 

The most significant phenotype of OE lines was upward 
curling of the leaf blades (Figure 4B). This phenotype ap- 
peared first at the seedling stage and then maintained 
through the rest of the plant growth (Figure 4B), The 
measurement of leaf rolling index (LRI) was performed at 
the time of flowering. The LRIs of the flag leaf, the second 



leaf from top, and the third leaf from top were approxi- 
mately 57%, 52%, and 50%, respectively (Figure 4C). The 
leaf blade cross-sections did not show significant differ- 
ences between wild type (WT) and OE in terms of 
organization of the sclerenchyma and vascular tissues. It 
suggests that there is no obvious alteration of leaf polarity 
in transgenic plants. However, in contrast to WT, the 
leaves of transgenic plants displayed smaller bulliform 
cells in the rolled regions (Figure 4D). 

Overexpression of 05MYBW3L affects the expression levels 
of several genes encoding cellulose synthase (CESA) 

To explore the possible molecular mechanism of 
OsMYBlOSL s effects on rice leaf development, we 
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Figure 4 Phenotypes and molecular analysis of OsMYBW3L overexpressing plants (OE). (A) OsMYBlOSL expression in OE as determined by 
qRT-PCR. Rice ACTINl gene was used as internal control. Each column represents an average of three replicates and bars indicate the SDs. 
(B) Morphology of 2-month-old WT and OE leaves. OE showed inward rolled leaves compared to WT (top panel); Cross sections of leaves of WT 
and OE (bottom panel). Bars= 10 cm (top panel) or 3 mm (bottom panel). (C) Statistical analysis of LRIs at the time of flowering for three OE lines 
of the top three leaf blades (n = 20). Asterisks {***) indicate a significant difference between transgenic plants and WT controls at P < 0.001, by 
Student's t-test. (D) Transverse sections of the flag leaves in the middle part of WT and OE. Shown are regions between two large veins. Iv, large 
vein; sv, small vein. Bars = 50 |jm. 



examined the expression profile in OE-1 and WT leaf 
blades using Digital Gene Expression profiling analysis 
(DGE). In all, 3,036 annotated genes showed expres- 
sions with more than twofold difference between OE-1 
and WT, of which 2,026 genes were up-regulated and 
1,010 genes were down-regulated (Additional file 1: 
Figure S2 and Additional file 2: Table SI). To obtain an 
ontological profile of the up- or down-regulated genes, 
ontological terms were assigned, and the enrichment 
significance was analyzed by AgriGO [51]. Twenty-one 
Gene Ontology (GO) terms were significantly enriched 
in the up-regulated gene set, showing that up-regulated 
genes were involved in multiple biological processes, 
such as metabolism, biogenesis, localization, and cellular 
processes (Figure 5A and Additional file 1: Figure S3). 
Meanwhile, we noticed that two GO terms, cellulose 
metabolic process (GO:0030243) and cellulose biosyn- 
thetic process (GO:0030244), were possibly related to the 
leaf rolling trait as reported in rll4 mutant [2] (Additional 
file 1: Figure S3). Therefore, we compared the expression 
levels of the rice CESA genes [52]. Our DGE data showed 
that OsCESAh OsCESA4, OsCESAV, OsCESAS, OsCESA9, 
and OsCESAll were up-regulated while OsCESAS and 
OsCESA6 were down-regulated (Additional file 1: Table S2). 
We verified the expressions by qRT-PCR for these 



CESA genes in three independent lines (OE-1, OE-4 
and OE-8). Four CESA genes, OsCESA4, OsCESA7, 
OsCESAS y and OsCESA9, showed significantly increased 
transcript levels; while one gene, OsCESA6, was down- 
regulated in the OE plants (Figure 5B). The qRT-PCR 
analysis was largely consistent with the DGE analysis 
except three genes, OsCESAl and OsCESAll of which, 
due to their relatively low Log2Ratio in DGE data. CESAs 
are proteins responsible for cellulose synthesis in plants 
[52-54], thus alteration of CESA expression may change 
the cellulose content. Therefore, we measured cellulose 
content in leaf blades of WT and OE plants. The results 
showed that OE leaf blades had 13% higher cellulose con- 
tent compared to WT (Figure 5C). These results suggest 
that OsMYBlOSL regulates CESA gene expression and 
may consequently influence cellulose synthesis in rice. 

Knockdown of OsMYBlOSL decreases the mechanical 
strength and cellulose content of leaves 

To explore the intrinsic function of OsMYBlOSL, we 
analyzed the phenotypes of mutants with decreased ex- 
pression levels. Since the loss-of-function mutants were 
not available, we created RNAi-knockdown plants for 
OsMYBlOSL. Six independent RNAi transgenic lines 
were obtained. Quantitative RT-PCR analyses showed 
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Figure 5 (See legend on next page.) 
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Figure 5 DGE, qRT-PCR and chemical analysis of wild type (WT) and OsMYBWSL overexpression (OE) plants. (A) Main clusters identified 
by GO enricliment analysis of genes up-regulated in OE compared with WT. In each GO biological process category, the bar indicates the fold 
enrichment, which is defined as -Log 10 (P value). The up-regulated genes are involved in multiple biological processes, such as metabolism, 
biogenesis, localization and transport. (B) Relative expression of rice CESA genes in WT and OE was determined by qRT-PCR analysis. The transcript 
levels of genes were normalized to the levels of ACTINl gene expression. The transcript levels of these genes in WT were arbitrarily set to 1. Values 
are the mean ± SD of three replicates. (C) Measurement of cellulose content in WT and OE. Cellulose content (milligrams per gram of total cell 
wall residues) of the flag leaf segments from WT and OE plants. The error bars were obtained from five measurements. Asterisks (**) indicate a 
significant difference between transgenic plants and WT controls at P < 0.01, by Student's t-test. 



that the expression level of OsMYBlOSL was significantly 
decreased in leaves of the RNAi lines compared with that 
of WT (Figure 6A). There were no obvious morphological 
changes in RNAi lines. However, in the heading stage of 
the RNAi plants, we noticed that the leaves were brittle. 
When the leaf blade was bent, WT leaves were flexible 
while the leaves of RNAi lines were broken easily, espe- 
cially in the midribs (Figure 6B). The force required to 
break the leaf blades was 50% less for the RNAi lines than 
that for WT (Figure 6C). The weakened mechanical 
strength may relate to the altered cellulose content in leaf 
blades. We measured the cellulose content in leaf blades 
of WT and RNAi plants. Compared with WT, the levels of 
cellulose content in RNAi lines were significantly de- 
creased (Figure 6D). Therefore, the decreased mechanical 
strength of RNAi plants might result from the decreased 
cellulose content. We determined the expression levels of 
ten CESA genes by qRT-PCR, and found that three CESA 
genes, OsCESA4, OsCESA7, and OsCESA9 were signifi- 
cantly down-regulated in the RNAi lines (Figure 6E). 
These results suggest that the down-regulated CESA 
genes in OsMYBlOSL RNAi plants may lead to a de- 
creased level of cellulose content and the weakened 
mechanical strength of transgenic leaf blades. 

Discussion 

By screening phenotypes of a transgenic rice population 
overexpressing rice genes in MYB family, we identified a 
rice leaf rolling-related R2R3-MYB gene, OsMYBlOSL, 
We showed its activity of transcriptional activation and 
described its influence on leaf rolling. We further showed 
that the altered expression of OsMYBlOSL resulted in 
changes in CESA gene expression, cellulose content, and 
mechanical strength in rice leaf. These results suggest that 
OsMYBlOSL may have important biological functions in 
rice leaf development. 

Leaf rolling is one of the frequently observed pheno- 
types in rice and other cereals. It is one of the important 
agronomic traits for high yield rice varieties [6-9]. Al- 
though many cultivars and mutants in rice show several 
types of leaf rolling phenotypes and a number of leaf 
rolling related genes have been identified [1-3,19,20,55], 
the mechanisms underlying the rolling phenomena re- 
main to be elucidated. It is generally realized that two 



major factors can influence leaf rolling. The first one is 
the establishment of ceU polarity and differentiation, 
and the second is physiological factors, such as osmotic 
pressure or turgidity in buUiform cells [1,2,56,57]. Previ- 
ous studies on mutants defective in leaf development in 
Arabidopsis thaliana, maize, and rice have shown that 
leaf rolling is related to the defective development of 
the leaf polarity along the adaxial-abaxial axis [1,58]. 
However, we did not observe obvious alterations in leaf 
polarity in the OsMYBlOSL overexpressing plants by 
histological analysis (Figure 4D) and SEM observations 
(Additional file 1: Figure S4). Instead, we found that the 
buUiform cells in the rolled regions of OE leaves were 
smaller compared with that of WT leaves (Figure 4D). 
Several investigations have shown that buUiform cells 
are involved in leaf roUing [2,59]. Loss of turgor pres- 
sure in buUiform cells could lead to leaf "roU up" during 
water stress [56]. Here the observed abnormal buUiform 
cells in the leaves of OsMYBlOSL overexpressing plants 
are consistent with those described in a number of known 
rolled leaf mutants [2,20] . For instance, cowl and rll4 mu- 
tants showed rolled leaves with serious shrinkage and 
smaller buUiform cells, but with no obvious defects in leaf 
polarity [2,20]. Therefore, our results suggest that the ab- 
normal buUiform cells, rather than defects in leaf polarity, 
may account for the rolled leaf phenotype in OsMYBlOSL 
overexpressing plants to some extent. 

Secondary cell walls constitute skeletal frameworks 
and furnish the plant body with mechanical strength 
[60,61]. Molecular and genetic studies in the model 
plant, Arabidopsis thaliana, revealed that a group of sec- 
ondary wall-associated NAC domain (SWN) transcription 
factors are master switches regulating a cascade of down- 
stream transcription factors, such as MYB transcription 
factors that leads to the activation of the secondary wall 
biosynthetic program [62,63]. It has been reported that 
AtMYBlOS {Atlg6S910) is directly activated by the mas- 
ter switches, SWNs [62,63]. In Arabidopsis thaliana, 
overexpression of R2RS-MYB gene AtMYBlOS led to a 
noticeable increase in the secondary wall thickness, while 
its dominant repression caused reduction in secondary 
wall thickness of both xylary fibers and interfascicular 
fibers [50,64]. It has been proposed that AtMYB103 has 
the potential to induce the expression of genes involved 
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Figure 6 Down-regulation of OsMYBIOSL decreases mechanical strength and affected cellulose biosynthesis. (A) Quantitative RT-PCR 
analysis for RNAi transgenic plants. The transcript levels of OsMYBIOSL in WT were arbitrarily set to 1 . Values are the mean ± SD of three replicates. 
(B) An easily broken OsMYBIOSL RNAi-1 leaves, as indicated by the arrow. (C) The force required to break WT and OsMYBIOSL RNAi lines flag 
leaves. Error bar represents SD (n =10 in each measurement). (D) Comparison of cellulose content between the wild type and OsMYBIOSL RNAi 
lines flag leaves. Error bar represents SD (n = 5 in each measurement). (E) Relative expression of rice cellulose synthase (CESA) genes in wild type 
(WT) and OsMYBIOSL RNAi plants was determined by qRT-PCR analysis. The transcript levels of genes were normalized to the levels of ACTINl 
gene expression. The transcript levels of these genes in wild type were arbitrarily set to 1. Values are the mean ±SD of three replicates. Asterisks 
or ***) indicate a significant difference between transgenic plants and WT controls at P < 0.01 or 0.001, respectively, by Student's t-test. 



in cellulose biosynthesis [50]. The transcriptional net- 
work regulating secondary cell wall biosynthesis may be 
evolutionarily conserved in plants, however, little is 
known about the detailed regulatory programs in grass 
species, such as rice [65]. 



Based on our analyses on the OsMYBIOSL and its 
encoded protein, we deduced that OsMYBIOSL might be 
an ortholog of AtMYBlOS in rice. We wondered whether 
it has the conserved functions in rice, a model plant in 
monocotyledon, and addressed this question by reverse 
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genetics, including overexpression and RNAi-knockdown 
approaches. Our study showed that overexpression of 
OsMYBlOSL resulted in rolled leaves and knockdown of 
OsMYBlOSL led to brittle leaves (Figures 4B and 6B). 
Further analyses disclosed that the expression of four 
rice CESA genes, OsCESA4, OsCESA7, OsCESAS, and 
OsCESA9, was increased significantly in transgenic 
overexpression rice plants (Figure 5B and Additional 
file 2: Table SI) whereas three of these CESA genes, 
OsCESA4, OsCESAV, and OsCESA9, were significantly 
down-regulated in the RNAi lines (Figure 6E). OsMYBlOSL 
is a putative R2R3 MYB transcription factor, which is 
expected to regulate gene expression by binding to the 
promoters of its downstream genes. Putative MYB bind- 
ing sequences (MBSs) can be identified by using the 
tools described in PLACE (A Database of Plant Cis-acting 
Regulatory DNA Elements; http://www.dna.affrc.go.jp/ 
PLACE/) [66]. We did find two and even more CNGTTR 
sites, the putative MBSs, in the respective promoter re- 
gions of these CESA genes (Additional file 1: Figure S5). 
Therefore, we hypothesized that these CESA genes would 
be the targets of OsMYBlOSL. Further investigation, such 
as chromatin immunoprecipitation (ChIP) and electro- 
phoretic mobility shift assay (EMSA), could be carried out 
to prove it. 

In higher plants, CESAs were characterized as catalytic 
subunits of cellulose synthase complexes [67,68]. In this 
study, our results showed that the cellulose content of 
mature leaves was higher in overexpression lines than 
that in WT (Figure 5C) while the OsMYBlOSL RNAi 
lines had decreased cellulose content in the leaf blades 
(Figure 6D). The results indicate that OsMYBlOSL may 
regulate the expression levels of CESA genes, which 
could subsequently alter the cellulose content. How- 
ever, we also noticed that the contents of cellulose were 
not positively correlated with the expression levels of 
OsMYBlOSL in transgenic plants, especially in the 
overexpression lines. Plant CESA genes are members of 
multigene families [69]. Arabidopsis has 10 CESA iso- 
forms which can be classified into six orthologous 
groups having non-redundant functions in cellulose 
synthesis as revealed by mutational analysis [5S] . Rice also 
has 10 CESA isoforms, of which OsCESA4, OsCESA? and 
OsCESA9 were showed non-redundant functions in cellu- 
lose synthesis specifically in secondary walls [52,70]. Based 
on these, we speculate that most rice CESAs might also 
play distinctive roles in cellulose synthesis. Some particu- 
lar CESA gene expression induction may not lead to sig- 
nificant change of cellulose content. Therefore, changed 
expressions of several CESAs in OsMYBlOSL transgenic 
lines may not induce obviously changes in cellulose con- 
tent of the transgenic plants. 

In addition, our study demonstrated that OsMYBlOSL 
RNAi-knockdown plants possessed brittle leaves with 



obviously decreased mechanical strength and reduced 
cellulose content in the leaf blades (Figure 6B-6D). Of 
note, reduced cellulose content and inferior mechanical 
strength are the common phenotypes of CESA mutants, 
such as irxl to irxS in Arabidopsis and brittle culmlS 
{bclS) in rice [70-72] . Thus, we wondered whether mech- 
anical strength was elevated due to the increased cellulose 
in OsMYBlOSL overexpression plants. Therefore, we also 
measured the mechanical strength of the overexpression 
lines of OsMYBlOSL, Nevertheless, the leaves of overex- 
pression plants did not show statistically significant 
changes in mechanical strength compared to wild type 
(Additional file 1: Figure S6). Mechanical strength is a 
complex trait, which can be affected by many factors 
[7S,74]. Perhaps, in our case, the increase of cellulose 
contents was not enough to cause the change of mech- 
anical strength in OsMYBlOSL overexpression plants. 

Moderate leaf rolling has been regarded as one of the 
important agronomic traits in high yield rice cultivars 
[6-9], and selection of this trait has been applied for rice 
breeding directly or indirectly in China [5,6]. One of the 
successful examples is the breeding of super-hybrid rice, 
Liang-You-Pei-Jiu (LYP9) [75,76]. LYP9 has moderate 
xrolling leaves and is widely cultivated throughout 
southern China due to its high yield and seed quality. 
The successful cultivation of LYP9 is therefore taken as 
an example for improving crop yield and quality using 
the rolled leaf trait [77]. Therefore, identification of 
genes that regulate leaf rolling would be beneficial to 
breeding. In this study, we report that overexpression of 
OsMYBlOSL in rice results in leaf rolling, which suggests 
that OsMYBlOSL could be a potential target for genetic 
manipulation of leaf shape and plant architecture in rice. 
In addition, leaf brittleness is an important agronomic 
trait for rice and other crops since it affects not only 
grain production but also the utilization of cereal straw 
as animal forage [74]. Here, we found that down- 
regulation of the expression level of OsMYBlOSL de- 
creases mechanical strength in rice plants, suggesting 
that OsMYBlOSL could be engineered for favorable leaf 
brittleness in rice, as well as in other cereal and forage 
crops. 

Conclusions 

In this study, an R2RS MYB transcription factor, 
OsMYBlOSL, has been identified in rice. We show that 
overexpression of OsMYBlOSL in rice results in leaf 
rolling, increased cellulose content, and elevated ex- 
pression of genes encoding cellulose synthase (CESA). 
Besides, the RNA interference (RNAi) -knockdown of 
OsMYBlOSL leads to a weakened mechanical strength 
and decreased cellulose content in leaves. Moreover, 
the expression levels of several CESA genes are up- 
regulated in OsMYBlOSL overexpression lines while 
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down-regulated in OsMYBl OSL-knockdown lines. Collect- 
ively, this work demonstrates that OsMYBlOSL impacts 
leaf shape, cellulose synthesis, and mechanical strength 
in rice. 

Methods 

Plant materials and growth conditions 

Rice {Oryza sativa L.) subspecies indica cultivar Kasalath 
and japonica cultivar Nipponbare were used for analysis. 
Rice seeds from control wild type and transgenic lines 
were immersed in water for 2 days and grown in a soil 
seed bed for approximately 30 days, and then the seedlings 
were transplanted to the field and/or greenhouse for 
further analysis. Plants maintained in greenhouses 
were grown under standard conditions at 28°C with 
16 h of light. The fields were located in the Experimental 
Stations (Beijing and Hainan) of the Institute of Genetics 
and Developmental Biology, Chinese Academy of Sciences, 
Beijing, China. 

Bioinformatics analysis 

A search for OsMYBlOSL homologs was performed 
using the National Center for Biotechnology Informa- 
tion (NCBI) BLAST server (http://blast.ncbi.nlm.nih.gov/ 
Blast. cgi). Domain prediction was performed using the 
Pfam database (http://pfam.sanger.ac.uk) and the NCBI 
Conserved Domains database (http://www.ncbi.nlm.nih. 
gov/Structure/cdd/cdd.shtml). Alignment was performed 
with ClustalW using default setting (http://clustalw.ddbj. 
nig.ac.jp/). Unrooted neighbor-joining phylogenetic tree 
was constructed using the neighbor-joining algorithm 
with PHYLIP version 3.572c [78]. Bootstrap analysis was 
carried out with 1,000 replicates. TreeView was used to 
generate the graphical output [79] . 

Subcellular localization 

For detection of subcellular localization of OsMYBlOSL, 
the open reading frame (ORF) of OsMYBlOSL was amp- 
lified from the FL-cDNA clone (AKl 11808) and then 
was cloned into vector pEZR(K)-LN to create the 35S: 
OsMYB103L-GFP vector. The 35S: OsMYB103-GFP and 
empty vector (pEZR(K)-LN) were transformed into rice 
protoplasts by polyethylene glycol-mediated transform- 
ation method [80] or onion epidermis by particle bom- 
bardment using a PDS-lOOO/He biolistic particle delivery 
system (Bio-Rad, California, USA), respectively. Subcel- 
lular distribution of the OsMYB103L protein was visual- 
ized by fluorescent confocal microscopy (Leica TCS SP5, 
California, USA). 

Transactivation assay in yeast cells 

The yeast strain Y2HGold, containing the AURl-C and 
MELl reporter genes with GAL4 binding elements in 
the promoters, respectively, was used as an assay system 



(Clontech, USA). The expression of AURl-C confers 
strong resistance to the highly toxic drug Aureobasidin 
A (Aba). The activities of a-galactosidase which MEL-1 
encoded were examined by X-a-gal staining. The full- 
length, N- terminal MYB DNA-binding domain (1-160 
Amino acids) or the C-terminal putative activation do- 
main (161-359 Amino acids) were respectively cloned 
into pBD-GAL4 vector. The pBD vector was used as a 
negative control. All these plasmids were individually 
introduced into cells of yeast strain Y2HGold. The yeast 
transformants were grown on SD medium in the ab- 
sence of tryptophan (SD/-Trp) for 2-3 days at 30°C. 
Transferred yeast cells were also grown on SD/-Trp in 
the presence of Aba and X-a-gal for 2-3 days at 30°C to 
monitor the generation of blue color. 

Promoter-GUS analysis 

The promoter of the OsMYBlOSL (about 1.4 kb upstream 
of ATG) was amplified from the rice genomic DNA (cv. 
Nipponbare) using the primers 5'-GGATCCACTCCTGC 
GAGGCTCTGATACG-3' and 5'-CCATGGATGCACT 
CT-AGATATCACTG-3'. Then, the DNA fragment was 
cloned into the pCAMBIA1301 (GAMBIA, Australia) 
vector resulting in fusion of the promoter and the GUS 
reporter gene. The ProOsMYBlOS: GUS construct was 
transformed into wfld type rice cultivar Kasalath. About 
20 independent transgenic lines were obtained after 
screening. GUS staining was performed according to 
the method of Jefferson et al [81]. 

Binary vectors construction and rice transformation 

For the overexpression construct of OsMYBlOSL, gene- 
specific primers, 5 -GGArCCTGCTAGCAGCTAGAT 
CAAG-3' and 5'-ACrAGrGTCATCCTCCTG-TGTT 
TATT-3' {BamH I and Spe I sites italicized) were used 
to amplify the coding sequence of OsMYBlOSL using the 
FL-cDNA clone (AKl 11 808) as template. The FL-cDNA 
clone was from the National Institute of Agrobiological 
Sciences, Japan [78] . The sequencing-confirmed PGR frag- 
ment was cloned into the vector pTCK303 [82] to create 
the overexpression vector of OsMYBlOSL driven by a 
maize ubiquitin promoter. 

For the RNA interference (RNAi) construct of OsMYBlOSL, 
a fragment of approximately 325 bp was amplified from 
OsMYBlOSL with two primers, 5' -GGTACCACTAGTCCG 
TGAAGCTGGCGATGAAC-3 {Kpn I and Spe I sites 
italicized) and 5'-GGArCCGAGCrCCACGACGAGCCC 
GAACAAC-3' {BamH I and Sac I sites italicized). The 
sequencing-confirmed PGR fragment was cloned into the 
vector pTCK303 as previously described [82] to create the 
RNAi vector of OsMYBlOSL driven by a maize ubiquitin 
promoter. 

These constructs were transferred to Agrobacterium 
tumefaciens strain LBA4404 by electroporation, and were 
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transformed into an indica rice cv. Kasalath, respectively, 
according to the method of Hiei et al [83]. 

RNA extraction and qRT-PCR analysis 

Total RNA was extracted from different tissues using 
Trizol reagent (Invitrogen, California, USA) according to 
the manufacturer s protocol After treatment with DNase 
to remove genomic DNA contamination, the first strand 
of cDNA was synthesized by M-MLV reverse transcript- 
ase (Promega, USA). Relative quantification of gene ex- 
pression by real-time PGR was performed on a Bio-Rad 
Chromo 4 Real Time PGR System (Bio-Rad, USA) with 
TransStart SYBR Green qPGR SuperMix Kit (TransGen, 
Beijing, Ghina). The housekeeping gene ACTINl was 
used as an internal control for normalization of RNA 
samples. Three replicates were carried out for each gene 
and each analysis was biologically repeated at least twice. 
Gene-specific primers are shown in Additional file 1: 
Table S3 line. 

Measurement of the leaf rolling index (LRI) 

To determine the leaf rolling index (LRI), two mea- 
surements were taken, Lw (Expand the leaf blade and 
determine the greatest width of the leaf blade) and Ln 
(Measure the natural distance of the leaf blade margins 
at the same site). LRI was calculated as LRI (%) = (Lw-Ln)/ 
Lw X 100 [3]. Data were collected on the flag (uppermost), 
second and third leaves on twenty individual plants at the 
time of flowering. 

Cellulose measurement 

Mature leaves were cut into pieces and ground into 
powder under liquid nitrogen. The materials were then 
extracted twice with 70% ethanol at 70°G for 1 h. Gell 
wall materials were dried under vacuum for cellulose 
content measurement [84]. Whatman 3MM filter paper 
was used to establish a standard curve for quantification 
of cellulose. The anthrone reagent was used to deter- 
mine the cellulose content. 

Breaking force test 

To examine the mechanical strength in leaf blade, the 
leaf blades of flag leaves from developmentally matched 
transgenic rice plants and wild type plants were collected 
and immediately used for measurement. The forces 
that break the samples at leaf blades were measured 
with a digital force/length tester (5848 Microtester; 
Instron, USA). 

Scanning electron microscopy (SEM) observation 

Samples were fixed in 2.5% glutaraldehyde solution. Fixed 
samples were dehydrated with gradual ethanol series, 
dried by critical-point drying method using liquid carbon 
dioxide (Model HGP-2, Hitachi, Tokyo), gold-coated with 



an Edwards E-lOlO ion sputter coater (Hitachi, Tokyo), 
and then observed using a S-3000 N variable pressure 
scanning electron microscope (Hitachi, Tokyo). 

Microscopy 

For light microscopic analysis, the samples were fixed in 
2.5% glutaraldehyde solution, and then dehydrated through 
a graded acetone series. The samples were embedded in 
Epon812 (SPI GHEM) and polymerized at 70°G. Sections 
were cut with a microtome (Leica RM2265, Germany) and 
stained with 0.25% toluidine blue, and finally observed 
with a microscope (OLYMPUS BX41, Japan). 

Digital gene expression profiling (DGE) analysis 

For digital gene expression analysis, transgenic overex- 
pressor (OE-1) and wild type flag leaves were harvested 
for RNA extraction using TRIzol (Invitrogen) as de- 
scribed by the manufacturer. The digital expression 
analysis was performed by the Beijing Genomics Insti- 
tute (http://www.genomics.cn) using the standardized 
procedure. The main reagents and supplies are Illumina 
Gene Expression Sample Prep Kit and Illumina Sequencing 
Ghip (flowcell), and the main instruments are Illumina 
Gluster Station and Illumina HiSeqTM 2000 System. All 
clean tags were mapped to the reference sequences (ftp:// 
ftp.plantbiology.msu.edu/pub/data/Eukaryotic_Projects/ 
o_sativa/annotation_Dbs/pseudomolecules/version_6. 1 / 
all.dir/all.cDNA) and only 1 bp mismatch is considered. 
Glean tags mapped to reference sequences from multiple 
genes were filtered. Remainder clean tags were designed as 
unambiguous clean tags. The number of unambiguous 
clean tags for each gene was calculated and then normal- 
ized to TPM (number of transcripts per million clean tags). 

The original data set is deposited in the National Insti- 
tutes of Health Gene Expression Omnibus database under 
accession number GSE52394. 

Accession numbers 

Sequence data used in this manuscript can be found in 
the Arabidopsis Information Resource (TAIR, http:// 
www.arabidopsis.org) and the rice genome annotation 
database (http://rice.plantbiology.msu.edu) under the 
following accession numbers: AtMYBlOS {Atlg63910), 
AtMYB26 {At3gl3890\ AtMYBSO {AtlgS7S60\ AtMYBSS 
{At4g01680), AtMYB61 {Atlg09S40), AtMYB67 {At3gl2 
720), AtMYB86 {At5g26660), OsMYB103L {OsOSgOSS 
20), Os01gS1260, Os01gl8240, OsOSg04820, Os07g314 
70, Os01gS0720, Os04g38740, OsActinl {Os03gS0890), 
OsCESAl {OsOSg08370), OsCESA2 {Os03gS9340), OsCE 
SA3 {Os07g24190), OsCESA4 {Os01gS4620), OsCESAS 
{Os03g62090), OsCESA6 {Os07gl48S0), OsCESA7 {OslO 
g32980), OsCESA8 {Os07gl0770), OsCESA9 {Os09g2S 
490), OsCESAll {Os06g39970), 
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Additional file 1: Figure SI. OsMYB103L-GFP is located to nucleus in 

onion epidermal cells. Figure 52. Number of differentially expressed 
genes in Digital Gene Expression profiling analysis (DGE) between wild 
type and OsMYBlOSL overexpressing plants. Figure S3. GO enrichment 
analysis of genes up-regulated in wild type and OsMYBlOSL overexpressing 
plants. Figure S4. Scanning electron micrographs (SEM) analysis. 
Figure S5. Schematic diagrams of the promoter regions of OsCESA 
genes. Figure S6. The force required to break flag leaves in wild type 
and OsMYBlOSL overexpression plants. Table S2. Expression levels of 
CESA genes in DGE analysis between wild type (V\^) and OsMYBlOSL 
overexpressing plants (OE-1). Table S3. Primers of qRT-PCR used in this 
article. 

Additional file 2: Table SI. Genes that up- or down-regulated two-fold 
in DGE data of OsMYBlOSL overexpressing plants compared to that of the 
wild type (Excel File). 
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